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CONDUCTION INVESTIGATIONS INTO THE MAGNETIC PROPERTIES OF MAHER 


By 

H. Chagini 1 and Robert L. Ash 2 
SUMMARY 

Because of problems resulting from poor quality control of the magnetic 
samples required for the original study, the graduate research associate has 
been able to study a different problem. This report presents the thermal 
analysis of a heat exchanger for heating air to temperatures on the order of 
3000 # C for use in NASA's high temperature wind tunnel. It was ascertained 
that an externally finned shell-and-tube type of heat exchanger with 
counterflow could be considered for this application. The methods of 
estimating the convective heat transfer coefficients are outlined in this 
report. As a part of thermal design the required size of heat exchanger was 
predicted. As a result of extreme length of heat exchanger ( = 1/2 kin), it 
was determined that a conventional heat exchanger may not be the most 
suitable tool for this application. 


Graduate Research Associate, Department of Mechanical Engineering and 
Mechanics, Old Dominion University, Norfolk, Virgin! a 23508. 

2 Chairman/Eminent Professor, Department of Mechanical Engineering and 
Mechanics, Old Dominion University, Norfolk, Virginia 23508. 



INTRODUCTION 

A neat excnanger is any device which provides the transfer of tnermal 
energy from one fluid to another. In the simplest type of heat exchangers 
the hot and cold fluids mix directly, as compared to more common types in 
which the fluids are separated by a wall. This type, called a recuperator, 
may vary in design from a simple plane wall between two flowing fluids to 
more complex configurations involving multiple passes, fins, or baffles 
(Ref. 1, 2). In this case conductive and convective and sometimes radiative 
modes of heat transfer must be considered. 

Many factors including thermal analysis, size, weight, structural 
strength, pressure drop, and cost, contribute to the design of a heat ex- 
changer. The heat transfer coefficient of each fluid is estimated oy tne 
geometry of the flow passages, fluid flow rates, temperatures and fluid 
properties. From these specified cond-itions an overall heat transfer co- 
efficient is then calculated. 

For our application it was determined to use an externally-finned shell 
and tube arrangement with counterflow to investigate the practicality of the 
task. Thermal analysis of this type of heat exchanger is presented in 
detail in the next section. 


THERMAL ANALYSIS 

The primary objective in thermal design of a heat exchanger is to de- 
termine the necessary surface area, A, required to transfer heat at a 
given rate for a given fluid temperature and flow rate. 

One of the first tasks in the thermal analysis of a heat exchanger is 
to evaluate the total tnermal resistance of the system. This is facilitated 
by employing the overall heat-transfer coefficient, U, defined as tne 



reciprocal of the sun of the thermal resistances of the system. For a 
cylindrical wall it can be written as (Ref. 3) 

U= (1) 

1/h. + [r. ln(r Q /r i ) ] /K +( -1) /h Q 

r o 

here h is the film coefficient, K is the tube thermal conductivity, r 
is the radius of the tube, and the subscripts o and i represent tne 
outside and the inside, respectively. 

To determine the rate of heat transfer, Q, from the energy balance 
between the two fluids 


Q = m C AT). . = m C AT) . . (2) 

M p 'hot p 'cold 

where m is the mass flow rate of flow, C is the specific heat at con- 

P 

stant pressure, and aT is the absolute value of the inlet and outlet 
temperature difference. 

Heat exchangers are designed by the basic equation (Ref. 4) 

Q = U A Tr 

For more complex type of heat exchangers, such as those involving 
multiple tubes, several shell passes or crossflow, the aDove equation is 
modified to include a correction factor (Ref. 4). 

The remaining term, AT, in the above equation is the logarithmic mean 
temperature difference (LMTD). The temperatures of the fluids in a heat 
exchanger are generally not constant, but vary from point to point as heat 
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flows. Even for a constant thermal resistance, the rate of heat flow varies 
along the path of the exchangers, since it depends on the temperature dif- 
ference oetween the hot and cold fluid. Figure 1 illustrates the change in 
tanperatures tnat may occur in a simple shell and tube counter flow heat 
exchanger. 



For a counter flow heat exchanger the IMTD is defined as (Ref. 3) 


AT a ‘ AT b 

In AT /AT. 
d D 


(4) 


based on the following assumptions: 

1. Rate of flow of each fluid is constant. 

2. Specific heat of each fluid is constant. 

3. Overall heat transfer coefficient through the exchanger is con- 
stant. 
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4. The system is adiabatic, heat exchange takes place only Detween the 
two fluids. 

5. Fluid flow is countercurrent. 

For the case of the finned heat exchanger, a convenient concept that 
can be used to provide a value for the heat transfer rate from a fin is the 
fin efficiency. The efficiency of the bar fin is given by (Ref. 3) 

4 .tanh_(pL)_ (5) 

PL 

Here p = h c /ky, L is the length of the fin and y is half the fin thick- 
ness. 

Use of Reynolds Analogy to Estimate 
the Convection Coefficients 

Empirical correlations for the heat transfer coefficients needed in 
this study were not available, out could be estimated using Reynolds ana- 
logy. By utilizing the boundary layer concept of no slip in the vicinity of 
a surface, the fluid is essentially stationary. As a result the transfer of 
heat takes place primarily by conduction. Reynolds analogy (Ref. 5) results 
from making the approximation that temperature and velocity cnange at pro- 
portional rates through the boundary layer as 

. q" u T 

- JL — / du = f dT (6) 

C t 0 'T 

p s s 

Here q" and x represent heat transfer and shear stress, respectively. 

Tne subscript s denotes the surface values. 

For the case of turbulent flow in a tube, integration of the above 
equation with upper limit of u = V and T = T^, where V is mean 
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velocity and is the bulk temperature, gives 


or 





( 7 ) 



(3) 


Shear stress at the surface, r , can be evaluated from the force balance 
on a cylindrical control volume of length L and diameter 0, as 


= (AP)D 
AL 


(9) 


but from Darcy's equation (Ref. 2), 

AP = f pVf_ 

L D 2 

where f is Darcy's friction factor. Therefore, 


( 10 ) 


final ly 


t $ « f V 2 /8 


St y = -JL = f/8 
pc p v 


( 11 ) 

( 12 ) 


where St^ is the local Stanton number. The same form is also valid for 
averaged values 
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In order to estimate the internal convection coefficient h an approx- 


imation of the friction factor f is needed. From Schilichting (Ref. 6) 
the friction factor, in terms of Reynolds number for smooth pipe, is given 
as 


— « 2.0 log (Re • f ) - 0.8 (14) 

f 


for our Re^ = 10 7 a friction factor of f * 0.008 was obtained. For cal- 
culating the convection coefficient parallel to the long pipe, tne flat 
plate approximation can be used. Similarly, the flow outside the pipe is 
turbulent, and from Reynolds analogy integration of equation (5) 


h 


x_ 

oC V 
P 



(15) 


where is the average skin friction. For a flat plate, can be 
approximated as (Ref. 7) 


1/5 

C f = 0.072(Re L ) 


(16) 


with Re^ * 5 x 10 6 an average skin friction of s 0.0035 was resulted. 



Calculation of Thermodynamic and 
Transport Properties of Flow 

Hot fluid is the result of the complete combustion of methane in the 


combustion 

chamber with excess oxygen 

in the process. The product 

mixture of 

the following gases. 



Gas 

Mole Fraction 


C0 2 

13.3% 


h 2 o 

26.5 


Air 

49.7 


0 2 

10.5 


At high pressure (= 100 atm) and temperature (mean average temperature 
of approximately 1300°C) the properties of pure gases are obtained from 
references 8 and 9. 

The coefficient of viscosity of a gas mixture of n components may be 
calculated by Wilke's semi -imperi cal correlation (Ref. 10) as 


v mix 


n 

E 

i = l 



n 


E 

j=l 



(17) 


where y.. are the viscosities of the component gases, 


x . and x . are 
i J 


mole fractions. The coefficients are given as functions of the vis- 

cosities and molecular weights of species i and j as 


[1 + (u-j/uj) ) ] 2 

. n i/2 


(18) 


[8(1 + M-/Mj) ] 


Similarly, the thermal conductivity of a mixture and polyatomic gases may be 
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divided into two portions as (Ref. 10) 


k . = k' . + k". 

mix mix mix 


(19) 


Here, k^ represents the monatomic thermal conductivity of the mixture, 

whereas kjj^ accounts for the diffusional transport of internal energy. 

Since the combustion product has no monatomic component the term k‘ . is 

mix 

zero, and k" . is also calculated from equations (16) and (17) as y. is 
’mix i 

replaced by the thermal conductivity of the gas components, K . (Ref. 11). 

A computer program was written in BASIC to calculate the viscosity and ther- 
mal conductivity of tne mixture. A listing of the program appears on the 
next page. 

The assunption of ideal gas made it possible to calculate the density 
of the mixture gas with the molecular weight of 28.25. Also, Prandtl number 
and specific heat at constant pressure were obtained based on the mole frac- 
tion of the mixture. 


PRESENTATION OF RESULTS 

Calculations are performed for heat exchanger tubes of various diam- 
eters of 5, 10, and 20 inches (0.125, 0.250, and 0.500 meters, respective- 
ly). The condition of staying within the compressibility limit of Mach < 
0.3 was maintained tnroughout the analysis. With sonic velocity of approx- 
imately 660 m/sec, a maximum speed of 200 m/sec is allowed. For the present 
analysis, it was assumed that tube thicknesses were Doth 1/16" (0.016 
meters). For all cases the fin length and thickness were taken as 0.016 
meters and 0.014 meters, respectively. They were assumed to be one fin 
thickness apart from each other. As a result a fin efficiency of 40 percent 
was calculated. 
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For the inner tube tne inlet and outlet velocities are calculated based 
on changes in the fluid density as a function of temperature. Methods of 
obtaining the convective coefficients have been already outlined in this 
report. As a main purpose of this study the required length and pressure 
drop across a smooth pipe as well as finned pipe heat exchanger were calcu- 
1 ated . 

Results for different tube diameters are presented in Table 1. 

Table I 


Inner 

tube 

diam (m) 

m 

max 

(kg/sec) 

No. of 
tube 

required 

Inner 

fluid 

tube 
velos . 

Convective 
coef f . 
(kv/m^* K°) 

u 

(Wm2-K°) 

Smooth 

Elpe 

mm 

length 

(m) 

Ap 

( atm. ) 

length 

(ra) 

Ap 

(atm. ) 

V . 
out 

V in 

HI 

m 

0.125 

27.3 

ll* 

190. U 

30.6 



1.97>* 

2699 

197.0 

901 

65.8 






(VI 

VO 






0.250 

100.7 

1. 

i66.y 

26.8 

♦ 

♦ 

1.925 

138U. 

39.7 

1*55 

13.1 






-3 







0.500 

L3L.6 

1 

167.0 

26.9 



1.898 

702 

1 9.9 

229 

3.2 


IlOTL': Outer tube diameters are twice the Inner tube diameters. 


CONCLUSION AND RECOMMENDATION 

Calculation of neat exchanger area for smooth as well as finned tubes 
indicates that a conventional heat exchanger is not the most feasible 
approach to this problem. Even if the heat exchanger pipes are assembled in 
a compact form of a loop, the very high pressure drop across the heat ex- 
changer cannot be tolerated. For the shortest length of finned tube in tne 
case of 20-inch pipe, a head loss of at least 3.2 atmospheres was calcu- 
1 ated. 

As an alternative, a systan which utilizes tne energy released from 
fusion process is recommended. This process can be adopted by applying 
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electrical current to a bank of sealed 
The energy released to the surrounding 
air needed for use in the wind tunnel. 


metal containers to melt the metal, 
air can provide the high temperature 
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